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Abstract 

We compute the time-dependent coverage in the random sequential 
adsorption of aligned d-dimensional cubes in R d using time-series expan- 
sions. The seventh-order series in 2, 3 and 4 dimensions is resummed in 
order to predict the coverage at jamming. The result is in agreement with 
Monte-Carlo simulations. A simple argument, based on a property of the 
perturbative expansion valid at arbitrary orders, allows us to analytically 
derive some generalizations of the Palasti approximation. 
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1 Introduction 



Random sequential adsorption (RSA) is a model for the random adsorption 
of non-overlapping objects on a surface. RSA is a conceptually simple model 
for irreversible behaviour. It has a rich variety of forms which are explained, for 
example, in the review paper of J.W. Evans Here, we consider the random 
filling of R d by non-overlapping aligned d-dimensional cubes. This is a simple 
model, whose one-dimensional version has been solved since a long time ||. In 
comparison, the level of understanding of this model in higher dimensions is low. 
It has been studied through 2, 3 and 4 dimensional Monte-Carlo (MC) simulations 
and through a 2-d perturbative analysis ||. A density series expansion || can 
also be used to analyze this model. It is equivalent to a direct perturbative 
expansion of the coverage in the time variable. 

A striking particularity of the model concerns the behavior of the fraction of 
space, 6d, filled at jamming. 6 d is numerically found to be close to the jamming 
coverage in one dimension raised to the d-th power. The relation 9d = 8^, 
proposed long ago by Palasti ||, is, in fact, one of the motivations of some MC 
experiments with an increasing precision. Nowadays, they allow us to rule out this 
exact relation. Nevertheless, as an approximation it appears to have a puzzling 
accuracy, e.g. of the order of 0.5% in two dimensions || 0. 

The Palasti relation is one of our main interests. More generally, we undertake 
the reconstruction of the coverage 6^(T) in 2, 3 and 4-d using as basic input its 
perturbative expansion (PE) in the time variable T. 

In the next section, this expansion is derived up to seventh order by using di- 
agramatic methods already given in Ref. f||, ||. The complexity of this expansion, 
even at moderate orders (6th or 7th), necessitates a computer algorithm. We will 
only give the results in 2, 3 and 4- (because the full analytical result in d is too 
long) . 

Using standard approximation methods, which combine PE and asymptotic 
(T — > oo) expansions, we extend the analysis of 9d of Ref.|| to 3 and 4 dimen- 
sions. This analysis is done in section 2, where it is found that 9d ~ Q\ with 
some decreasing precision as d increases. In section 3, we provide a functional 
understanding of the Palasti relation as a consequence of a particular property of 
the PE, valid at arbitrary orders. Finally, by taking into account the asymptotic 
regime, we give some improved generalizations of the Palasti approximation. 

2 The d-dimensional perturbative expansion 

We will study the RSA of hypercubes in the continuum limit on R d . This 
is defined as the scaling limit of RSA of cubes on a d-dimensional hypercubical 
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lattice. The discrete model, in which the size of the adsorbed cubes enters as a 
parameter, has its own interest and it has been also investigated ]2j, [|, [7], |j in 
1 and 2-dimensions. At small sizes, the Palasti approximaton works perfectly, 
and more generally, the model is easier to handle [§, |10| . We consider here the 
discretisation as a convenient regular isat ion of the PE from which we will keep, 
order by order, the scaling limit. 

We follow the definition of the discrete model given in Ref. f7|. The adsorbed 
cubes have their edges parallel to the d axes of the lattice, and they are made of 
k d unit cells, k being an integer. Their volume l d = (kb) d enters in the definition 
of the dimensionless time, T, T = ptl d , in which b is the lattice spacing, p is 
the attempt rate to adsorb cubes per volume unit. On the lattice, the number of 
attempts per site after time t is thus ptb d , i.e. T/k d . Therefore the coverage can 
be written as 

e d (k,T) = k d Pd (k,^) (l) 

The coverage 9 d {T) in the continuum limit is then given by the limit of 8 d (k,T) 
as k goes to infinity with T fixed. 
The PE of p d (k,t) is: 

(-t) n 

p d (k,t)=t-J2Cn(k,d) (2) 

n>2 n - 

where the first order in time is t, because the first adsorption attempt on an 
empty lattice is always accepted. We refer to Ref.|| for a comprehensive and 
complete derivation of the PE. It starts with an operator formalism of the RSA 
mechanism and ends up with diagrammatic rules, that are sketched below, for 
computing the coefficients C n (k,d). 

To generate the n-th order diagrams, one expands the expression: 

na-na-^i)) 

i=2 j=l 

The set of variables defining a deposition is denoted by % , and Ki j is a hard 
Mayer function (Ki j = 1 only if i and j are overlapping objects, otherwise). All 
monomials of this polynomial are connected labeled graphs which are regrouped 
in classes of the same topology. Thus, they appear with combinatorial weigths 
w nt i in such a way that 

C n (k,d) =^2w n , ti I(T nA ) . (3) 

i 

I(^ n ,i) is th e graph contribution which embodies the details of the model. 

In order to compute I(T nti ), one has to do the summation over {xi}, the n 
vertices of T ni . At this stage, the vertices, {x;}, are to be understood as d- 
dimensional lattice points, i.e. xi = (xj,xf, . . . ,xf), because a d-dimensional 
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aligned hypercube can be characterized by the position of a distinguished point, 
for example a corner. Finally, the graph contribution is 

J(rn,i) = E U K ( X ^ X rn p )=In,i{k,d) . (4) 
{x} V 

The product flj is done on the links of Y n ^ and the sum is over all the lattice 
degrees of freedom for every variable xi, except one which can be frozen at the 
origin by transational invariance. 

Any graph contribution, I ni (k,d), is thus the d-th power of the contribution 
of the same graph in one dimension. 

I n 4k,d) = I n ,(k,l) d . (5) 

This simple property, that we shall use in the next section, is a consequence of the 
construction rule and of the factorization property of hard Mayer functions 
for hypercubes: K d (x,y) = Yu=i Ki(x\ y l ). 

Finally, in the continuum limit, the PE of 6d{T) is, from (JJ) and (0): 

6 d (T) = T - C n (d) T n (6) 

where, from (|H|) 

C n (d) = lim{A; 1 - n J n ,<(^l)} d = E^ l nt (7) 

The coefficients C n (d) are well defined, because I n! i(k, 1) are polynomials of degree 
(k — 1) in the variable k that we have analytically computed. Their leading 
coefficients are denoted by J n> j in Eq.(^). 

We have thus, stored in our computer algorithm, all the J n i and w n ^ necessary 
to know C n (d) analytically in d, up to the seventh order in T. Nevertheless we do 
not explicitly write them here due to their number (for n = 7 there are about 10 3 
graphs). At fixed d, C n (d) are just the rational numbers that we give in Table 
I for d — 2, 3 and 4. We have performed various tests of our expansion. The 
simplest is for d = where (T) = 1 — e~ T . For d — 1 the exact solution is 
known for any k and in the scaling limit 9\{T) reads: 

9 1 (T) = [ T dte~ 29{t) with g(t) = /* dx 1 ~ & * . (8) 
Jo Jo X 

The agreement was checked analytically in k, including the scaling limit. Finally, 
in 2-d, we recover the expansions given for k = 2 and k = oo in Ref.||. 

These coefficients C n (d) can be used in various ways to approximate 9d{T). 
We give, in the following, our analysis. In addition to the PE Eq.(|^), one knows 
how the jamming limit, 6d = 9d(T = oo), is approached [[□], 



9 d (T)-6 d = 0{\nT d - 1 /T] (9) 
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Standard methods which embody these two kinds of information have shown, in 
other contexts, their ability to reasonably approximate the coverage even in its 
asymptotic regime [Q, ||. 

We give, as a first attempt, the predictions for 9d from a staigthforward ex- 
tension of the method used in Ref.||. One uses the mapping a>(T), 

1 + In (1 + ^-T)^ 1 

"CO = i T+w ' b - 1 (10) 

which behaves as T for small T and as 1 — \nT d ~ 1 /bT for T — > oo, to convert the 
PE, Eq.(|]), into a power series of uj. The Pade approximants formed from this 
series, when evaluated at u — 1, give a ^ consistent with the behaviour (H), for 
any choice of the b which is used as a variational parameter. 

Let us first recall the 2-d results from Ref.||. Among the sixth and sev- 
enth order approximants, [5,2], [3,4], [4,2] and [3,3] cross near b = 1.30 yielding 
9 2 = 0.52625(5), in agreement with the data of Ref.g, 9 2 = 0.562009(4). As d 
increases, we expect some continuity, and we indeed observe a crossing at b = 1.40 
for the [4,2] and [3,4] approximants, and for the remaining ones [5,2] and [3,3] 
there is a crossing at b = 1.45. Nevertheless, the prediction, now, loses its accu- 
racy, because these intersections give 9 3 = 0.40 and 0.46, respectively. Moreover 
in 4-d this intersections disappear and the best determination of b is b ~ 2.3 — 2.5 
At this value, the [4,3], [4,2], [3,3] and [3,2] Pade-approximants give #4 = 0.30(1) 

These values are consistent with the MC data since #3 = 0.4227(6), 0.4262 
0.430(8) and 0.422(8) according to the Refs.|L2L 0, |l|, ||] respectively and 6 4 = 



0.3129 or 0.3341 from Ref.[||] and Ref.|T3| respectively. Obviously, the Palasti 
values, 9f = 0.7476, 0.5589, 0.4178, and 0.3124, d = 2,3,4 cannot be discussed 
(except in 2-d) within the poor precision of the previous results. Nevertheless, 
we show in the following section that perturbation theory is able to provide a 
functional insight into the Palasti approximation. It goes beyond a numerical 
coincidence. 



3 Variational approach to the Palasti approxi- 
mation 



The RSA mechanism of deposition of aligned cubes is characterized by power 
law behaviors in d for both the coverage at jamming, in the Palasti approximation, 
and the graph-integrals, in the framework of the perturbative expansion. We 
shall give in this section a variational derivation of this puzzling approximation 
together with an extension to finite values of T . 

Starting with the representation of d (T) given by Eqs.(||[7|), 

6 d (T) = £ A hn I hn d T n (11) 

i, n>l 
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we write, in order to emphasize the role of the power law behavior (|5|), that 

d-l j\ 

h% = ikn - ab n ) d + E(-1)^ +1 37T— ^ b n{d - r) V (12) 

r=0 \ r 

This is nothing but a binomial identity, valid for any a and b. We can insert (p^) 
into ([ll]) t° obtain a two-component representation of 9 d (T): 

6 d {T) = 6 d (T) + e d (T) (13) 

where 

d-l j| 

^) = E(- 1 ) d ' r+1 F73 ' „ ^' r f r (6^T) (14) 
^ r!(d-r)! 

e d (T) = E Ai, n (Ji, n - ab n ) d T n (15) 

i, n>l 

In the representation (|I3|-|T4|), 9 d (T) involves all the lower dimensionality cov- 
erages, and the a priori d- and T- dependent parameters a, b. A natural way 
to choose of a and b is thus to look for maxima of 6d{T), i.e. minima of the 
perturbative component e d (T). 

Beginning with the 2-d case, one finds that #2 CO = 2a@i{bT) — a 2 6 (b 2 T), has 
a maximum with respect to a : 

OUT) - <^ T > 



for a = 9i(bT) / 9o(b 2 T) , i.e. a value of a where e\ vanishes. If e% is , this is 
nothing but the Palasti approximation at T = 00. 

This result can be extended to any d: Under the assumption e r = for 
1 < r < d — 1, Od{T) has an extremum 0* d (T): 



e* d (T) 



efib^T) eUibT) 



Oq (b d T) 9 d _ 2 (b 2 T) 



for a = O^-^/e^T). 

This is simply because the minimizing condition can be reexpressed in terms 
of e^-i if the a parameter is assumed to be ^-independent. The relation is: 

de d( T ) . fur \ 

— = -d e d -i{bT) . 

oa 

We are thus lead to write the PE flT3|) of 9 d (T) under the form: 

6 d (T) = 6 d (T) + e d (T) (16) 

where we choose 

^ T > - (17) 
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as a convenient generalization of the Palasti approximation on the whole time- 
range. To go further, one can eliminate T, order by order, in terms of 9* d (T) 
(9* d (T) ~ T at T = and 9* d (T) increases if b > 1), and then, substitute this 
expression in €d(T), (ed(T) ~ T 2 ). One finally obtains 9d(T) as a power series 
in 9 d (T), i.e. as a perturbative expansion where the first order is the Palasti 
approximation. 

In this approach, it remains to take into account the asymptotic behaviour 
(0). One oberves that ^(T6^(T)) and 9 d (T) reach their asymptotic limit in the 
same way, 0[ln T d ~ 2 /T}. On the other hand, ^(T8 d (T)) is perturbatively of the 
form (0), so that the previous arguments apply. Thus we propose: 



as an approximation of Od(T) with a minimal component e^(T). The jamming 
coverage is then 9 d = 9 d (oo) + e d (oo) , where 



We have to verify, in PE, that €d{T) is small and that a precise determination 
of the b parameter can be done. We express ed(T) as a power series in 6* d {T) 
and form its Pade table. We find e^(T = oo) from these approximants when 
9* d = 8^(00). Here, b is used as a variational parameter. 




8* (00) = 81/8, 

ei(oo) = 8*J 82 



1 = 



0.4225 



0.3176. 
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In 2-d, 3-d and 4-d, we find a minimal dispersion of the [3,4], [4,3], [5,2], [4,2], 
[3,3], [2,4], [2,3] and [3,2] Pades at the values: 



d = 2 b ~ 1.32 e 2 (oo) = 0.004(2) 
d = 3 6=1.35(3) e 3 (oo) = 0.00(1) 
d = 4 6=1.40(5) e 4 (oo) = 0.01(1) . 

6d(T) is found to be positive on the whole physical region (except near T = oo for 
some marginal solutions in 3-d) and is of the order of 6^(00). We finally obtain 
that 

9 3 = 0.42(1) and 6 A ~ 0.32(1) . 

In most of the cases, the positivity of e^(T) indicates the violation of the Palasti 
conjecture. 



4 Conclusion 



This study has been devoted to the RSA coverage of R d by <i-dimensional 
aligned hypercubes, especially in 3 and 4 dimensions. For these dimensions, we 
have explicitely presented the perturbation result to seventh-order. The series 
in other dimensions and any discrete case are available upon request. The con- 
tinuum model at long times is expected to be more difficult to describe than its 
discrete analogue. This is supported by the fact that we cannot predict 6d be- 
yond a moderate accuracy with our perturbative information. Nevertheless, we 
indicate, in this framework, how one can approach the Palasti conjecture and its 
corrections. 

To conclude, we want to mention that the simple rule 6 d ~ @d-i/@d-2 applies to 
various discrete models, in particular RSA on rf-dimensional lattice with nearest- 
neighbour exclusion |16|], where #0 = 1/2 and 9\ = (1 — e~ 2 )/2. As d runs from 2 



to 6, we obtain the sequence of approximations (0.374, 0.307, 0.254, 0.229, 0.206) 
which are quite close to the corresponding data (0.364, 0.304, 0.264, 0.233, 0.209). 



7 



References 

[1] J.W. Evans, Review of Modern Physics (to be published). 

[2] P.J. Flory, J. Am. Chem. Soc. 61, 1518, (1939). 

J.J. Gonzalez, PC. Hemmer and J.S. Hoye, Chem. Phys. 3, 228, (1974). 
A. Renyi, Sel. Transl. Math. Stat. Prob. 4, 203, (1963). 

[3] R. Dickman, J.S. Wang and J. Jensen, J. Chem. Phys. 94(12), 8252, (1991). 

[4] G. Tarjus, P. Schaaf and J. Talbot, J. Stat. Phys. 63, 167, (1991). 

[5] I. Palasti, Publ. Math. Inst. Hung. Acad. Sci. 5, 353, (1960). 

[6] B.J. Brosilov, R.M. Ziff and R.D. Virgil, Phys. Rev. A 43, 631, (1991). 

[7] V. Privman, J.S. Wang and P. Nielaba, Phys. Rev. B 43, 3336, (1991). 

[8] J.W. Evans, Phys. Rev. Lett. 62, 2642 (1989). 

[9] R.S. Nord and J.W. Evans, J. Chem. Phys. 82 (6), 2795, (1985). 

[10] P. Schaaf, J. Talbot, H.M. Rabeony and H. Reiss, J. Chem. Phys. 92, 4826, 
(1988). 

[11] R.H. Swendsen, Phys. Rev. A 24, 504, (1981). 

[12] W.S. Jodrey and E.M. Tory, J. Stat. Comp. Sim. 10, 87, (1980). 

[13] B. Blaisdell and H. Solomon, J. Appl. Prob. 19, 382, (1982). 

[14] D.W. Cooper, J. Appl. Prob. 26, 664, (1988). 

[15] R.S. Nord, J. Stat. Comp. Sim. 39, 231, (1991). 

[16] A. Baram and D. Kutasov, J. Phys. A: Math. Gen. 22, (L251), 1989. 



8 



cm 


n 


d=2 


d=3 


d=4 


2 


4 


8 


16 


3 


23 


101 


431 


4 


168 


27 


A QOflQO 

27 


5 


105895 


30622397 


7785588283 


72 


864 


10368 


6 


6709687 


11820085427 


17179096575091 


450 


13500 


405000 


7 


1385692277 


5815177641671 


183314989218969359 


8100 


233280 


65610000 



Table I: Coefficients of the time-series expansion of the coverage 9 d (T) defined 

by 

(-l) n 

e d {T) = t-Y, K —r c n (d) T n . 

n>2 U] 
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